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scalar fields minimally coupled to the Einstein-Hilbert gravity and argue that the counting 
of degrees of freedom can be done for scalar theory and gravity separately. In this paper 
we investigate the system with only two Stiickelberg scalar fields. In this case we find the 
analytic expression for the determinant of the kinetic matrix of the scalar field Lagrangian 
and perform the full constraint analysis. In 1 + 1 space-time dimensions the theory cor- 
responds to the full non-linear massive gravity, and this determinant vanishes identically. 
In this case we find two first-class constraints, and present the corresponding gauge sym- 
metry of the theory which eliminates both scalar degrees of freedom. In 3 + 1 dimensions 
the determinant of the kinetic matrix does not vanish identically and, for generic initial 
conditions, both scalar fields are propagating. 
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1 Introduction 

The observation of the accelerated expansion of our universe is the driving motivation for 
various infrared modifications of general relativity. One of the theoretically most natural 
infrared modification would be to give a small mass to the graviton. Since the early discov- 
ery of the quadratic Fierz-Pauli mass term for metric perturbations in [1], there has been 
an ongoing search for a healthy non-linear completion of massive gravity. The construction 
of the non-linear graviton mass term is based on the use of an auxiliary non-dynamical 
reference metric, which as an absolute object would break the diffeomorphism invariance 
of general relativity. The diffeomorphism invariance can be restored by introducing four 
Stiickelberg scalars, corresponding to the four coordinate transformations [2-4]. However, 
a generic theory of four Stiickelberg scalars together with the two degrees of freedom of 
massless graviton propagates six degrees of freedom in total. It is one degree of freedom 
too much in comparison to the five degrees of freedom expected from the massive spin-2 
representations of the Poincare group. Moreover, the additional degree of freedom is sick 
and represents the (in)famous Boulware-Deser (BD) ghost [5]. 

After an order-by-order construction of a non-linear theory which is ghost-free in the 
decoupling limit in [6] , a full resummed theory of non- linear massive gravity was proposed 
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by de Rham, Gabadadze, and Tolley (dRGT) [7]. In unitary gauge this theory has been 
shown to propagate five degrees of freedom [8, 9]. The Hamiltonian analysis of the full 
diffeomorphism invariant theory including the four Stiickelberg fields also seems to confirm 
the expectation that the dRGT theory propagates at most five degrees of freedom [10-12] 
(for recent counterarguments see [13]). However, the canonical analysis of dRGT theory 
in the presence of the four scalar fields is intricate, and in the existent literature it is often 
obscured either by mixing the gravitational and scalar degrees of freedom or by introduction 
of new auxiliary fields. 

In the present paper we take a different point of view and treat dRGT massive gravity 
as a theory of Stiickelberg scalar fields (f) A coupled to the Einstein-Hilbert gravity. Since 
the theory is reparametrization invariant, and the scalars are coupled to gravity minimally, 
we shall count the degrees of freedom propagated by the metric and by the scalar fields 
separately. Hence the absence of the sixth mode in dRGT theory should manifest itself as 
the feature of the scalar fields Lagrangian alone. 

Motivated by these considerations we study the dynamics of the Stiickelberg scalar 
fields given by the dRGT mass term [7]. We observe that, if seen as a particular scalar field 
theory, the dRGT scalar field Lagrangian allows for an arbitrary number of scalar fields in 
it. In particular, the number of scalar fields N can be chosen to be less than the space-time 
dimension d+1 without affecting the diffeomorphism nor the space-time Lorentz invariance 
of the theory. We dub the dRGT theories of gravity with reduced number N < d + 1 of 
Stiickelberg scalar fields as "reduced massive gravity" . 

The simplest particular cases of such dRGT inspired scalar theories include, for d = 0, 
the action of a massive relativistic particle in N dimensions and, for N = 1, the single 
"k-essence" field with DBI-like action [14]. Another "simple" choice is arbitrary iV fields 
in 1 + 1 dimensions, and gives the action of a relativistic string in TV-dimensional target 
space-time. In the case iV = 3, with three scalar fields living in a configuration space 
diffeomorphic to M 3 , the reduced dRGT action can be regarded as a particular effective 
field theory of homogeneous solid [15]. The degree of symmetry of the solid depends on 
the isometries of the metric /ab(</>) in the internal space of scalar fields. If the metric is 
symmetric under the SO(3) group, and the action contains only the term, invariant under 
the volume preserving diffeomorphisms, then it describes a perfect fluid. The case with 
the number of scalar fields N > d + 1 has been recently discussed in [16, 17] as a theory of 
multiple Galileon fields covariantly coupled to the dRGT massive gravity. 

In general, the solutions of the reduced massive gravity theories are expected to break 
Lorentz and rotational symmetries and lead to anisotropic cosmologies. The pattern of such 
breaking is determined by the number of scalar fields and the signature and isometries of 
the reference metric. The connection of reduced massive gravity theories to the Lorentz 
violating massive gravity theories will be discussed in more detail in the main body of the 
paper. Another possible application of reduced massive gravity theories could be found in 
modeling the translational symmetry breaking and momentum dissipation in holography. 
In particular, in [18] the conductivity in the boundary theory was calculated in the presence 
of a Lorentz violating graviton mass term in the bulk, that originated from the dRGT-like 
action with two Stiickelberg fields and Euclidean reference metric. The models discussed 
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in our paper could be further used in holographic constructions. 

In this paper we consider the case of reduced massive gravity with two Stiickelberg 
fields. It is the simplest case with several scalar fields involved, in which we can write the 
Hamiltonian and constraint structure explicitly. We perform the full Hamiltonian analysis 
of the scalar field sector and find that, in distinction from the dRGT massive gravity the 
determinant of the kinetic matrix does not vanish. Hence the scalar field Lagrangian in 
general propagates two degrees of freedom. We formulate the condition for the scalar field 
configurations on which the determinant vanishes and investigate the different regions in the 
phase space of scalar fields. We show that on the singular surface, where the determinant 
of the kinetic matrix vanishes, the theory is equivalent to 1 + 1-dimensional massive gravity 
and thus has no dynamical degrees of freedom. We also show that the regular solutions 
away but in close vicinity of the singular surface approach the singular surface but can 
never reach it in finite time. At the same time any perturbation of the singular solution 
drives the system away from this singular surface. In quantum theory the vanishing of the 
determinant signals the strong coupling regime for the scalar fields, and the dynamics in 
the vicinity of the singular surface are highly affected by quantum corrections. Whether 
or not the two dynamical degrees of freedom away from the singular surface contain ghost 
modes might depend on the particular choice of the reference metric in the configuration 
space of the scalar fields. We do not address this question in the present paper, but leave 
it for future studies. 

The paper is organized as follows. In section 2 we recall the formulation of dRGT mas- 
sive gravity. In section 3 we formulate the theory of reduced massive gravity and perform 
the Hamiltonian analysis away from the singularity surface. In section 4 we consider the 
behaviour of the system on the singular surface, and show that it is equivalent to 1 + 1 
dimensional massive gravity. We perform the canonical analysis in this case and find the 
gauge symmetry of the scalar fields, eliminating both scalar degrees of freedom. Section 5 
is devoted to conclusions. 

2 Non-linear massive gravity in Stiickelberg formulation 

The non-linear massive gravity action can be written in terms of the variables 



where g^ p is the inverse space-time metric, and f^u is an auxiliary reference metric. The 
full dRGT action is given by 







71=0 



where the characteristic polynomials e n (X) of a 4 x 4 matrix X are 
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The squared brackets denote the traces, and the coefficients a n are arbitrary. It is also 
possible to rewrite the mass term in terms of the characteristic polynomials of the square 
root matrix (^jTCj = yy/g^fj as 

4 

= V=gJ2^^ . ( 2 - 3 ) 

k=0 

with the coefficients /3& given by 

*=B-»*(;:i)*.- <"> 

The characteristic polynomials of an n x n matrix X can be rewritten as the characteristic 
polynomials of its eigenvalues Aj as e„(X) = e n (Ai) [19, 20], where 



eo 



(Ai) = l, e 1 (A,) = ^A i = [X], 



e2(Ai) = ^ XiXj 

i<j 



e n (Ai) = A1A2 . . . A„ = detX . 



Since here the matrix X = yfi is a square root matrix, we note that its eigenvalues are, 
by definition, equal to the square root \/Ai of the eigenvalues of the matrix Q. Hence the 
mass term (2.3) can be rewritten in terms of the eigenvalues of the matrix Q„, without the 
need of finding the explicit expression of the square root matrix itself, as 



= v^X^^V^) • (2-5) 

fc=0 

Since the mass term (2.3) explicitly depends on the auxiliary metric fav it breaks the clif- 
feomorphism invariance of general relativity. It can be fully restored by introducing four 
Stuckelberg scalar fields <f> A , ^4 = 0, 1,2,3, corresponding to the four coordinate transfor- 
mations as fav = du</> d u <fi B rijiB [3]. In addition, in this parametrization the auxiliary 
metric is invariant under the Lorentz transformations in the scalar field space [4]. 
Hence the scalar field indices A, B are raised and lowered with the Minkowski metric 
rjAB = diag ( — + + + ). In this case the reference metric / M „ is said to be 'flat' since it 
is simply a coordinate transformation from the flat Minkowski metric t)ab- An arbitrary 
'curved' reference metric fav can be obtained by replacing the flat metric t]ab with some 
arbitrary scalar field metric / J 4s(</ ) ) [21]. Then the Lorentz transformations in the scalar 
field configuration space are replaced by the isometries of the metric /Us (</>)• 

The Stuckelberg formulation of the massive gravity allows for an equivalent form of 
the mass term (2.3) by introducing a diffeomorphism invariant matrix 

l£ = gTd^du^fBC ■ (2-6) 
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Since the traces and eigenvalues of the matrices Xg and Q„ = g fJ,p d p (f) A d l/ (f) B Jab are equal 
then the mass term (2.5) can be equivalently written in terms of the eigenvalues of X^. 
This rewriting makes manifest that any non-linear massive gravity theory can be viewed 
as a theory of a number of scalar fields minimally coupled to gravity. 

3 Reduced massive gravity 

In the present paper we adopt the point of view that the mass term Lagrangian used 
in the non-linear dRGT massive gravity is a Lagrangian describing four Stuckelberg scalar 
fields coupled to gravity. The motivation of restricting the number of scalar fields to the 
number of space-time dimension in the context of non-linear massive gravity is that around 
the background solution g pv = 77^, <f> A = x^8 A the metric perturbations have a Lorentz 
invariant mass term of the Fierz-Pauli form at the quadratic level. However, if seen as 
describing a theory of scalar fields, the action 

4 

= V^Y, T b= S^d^du^fBcW (3-1) 

n=0 

describes just some particular theory of derivatively coupled scalar fields, and depends only 
on their first derivatives. This theory is diffeomorphism invariant even when the number 
of scalar fields is not equal to the space-time dimension. Therefore from the scalar field 
theory point of view the number of scalar fields N can be chosen arbitrary, both less or 
greater than d+ 1. In the case N 7^ d+1 the matrices Xg and Q„ have different dimensions, 
N xN and (d+1) x (d+1) respectively. Nevertheless, the non- vanishing eigenvalues of these 
matrices are equal, and both formulations (2.2) and (3.1) of the action are still equivalent, 
even though the formulation in terms of the smaller matrix is evidently simpler. 1 

In this work we focus on the case of two scalar fields cj) A = {<f>°, (f) 1 } in 3 + 1 dimensions 
as the simplest non-trivial case inspired by the dRGT massive gravity. Using the diffeo- 
morphism invariant variables X^ proves to be particularly useful for this setup since the 
matrix X^ is a 2 x 2 matrix in this case whereas 0^ in a (3 + l)-dimensional space-time is 
a 4 x 4 matrix. The action of the scalar fields in any d > 1 then takes the simple form 

^ = («o + «iTr(I - VX) + a 2 det(I - VX")) , (3.2) 

where we have used the fact that for any 2x2 matrix X the polynomials 63,4 (X) vanish. 
In the case ao = ot\ = and for the scalar field metric taken to be the Minkowski metric 
rjAB, the full theory C = Ceh + m 2 £-<f> has the solution 

tf a, = rT, (f>° = x°, ^ = x\ (3.3) 
The quadratic action for the perturbations 

^ = X A=(^_ X ^ 5 A + ^X_ X ^ 5 A (34) 

1 The coefficients a n , a n in (2.2) and (3.1) respectively coincide only when the number of fields equals 
the space-time dimension, i.e. when N — d + 1. 
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then reads 

Cf = 2 [(Z* 01 ) 2 - h^h 11 ) + 4h AB [r, AB d c X C ~ d AX B] +2 [{d AX A f - d A X B d BXA ) , (3.5) 

and the indices A, B = 0,1. In 1 + 1 space-time dimensions this action coincides with 
the action for metric and scalar field perturbations around the Minkowski background in 
massive gravity. However, in 3 + 1 space-time dimensions, this corresponds to a Lorentz- 
violating Fierz-Pauli-type mass term for metric perturbations. A thorough analysis of 
Lorentz-violating graviton mass terms, preserving the Euclidean symmetry of the three- 
dimensional space was carried out in [22] (see also an earlier work [23]). In our case, the 
symmetry of the three-dimensional rotations is in general not preserved by the ground state 
of the theory. Therefore, the possible mass terms of our theory go beyond the considerations 
of [22]. Out of their investigated mass terms, only the mass term with ni2 = ms ^ 
0, Tio,i,4 = can be obtained in reduced massive gravity (if the number of scalars N = 3). 
However, the stability analysis of [22] does not directly apply to our case since the number 
of Goldstone fields is different. 

A particular instance when the (3+l)-dimensional dRGT theory of massive gravity 
reduces to the special case of two Stiickelberg fields is the case of a degenerate refer- 
ence metric. To see this one can consider the spherically symmetric ansatz <fp = f(t,r), 
4> % = g(t, r), i = 1, 2, 3. 2 For the flat auxiliary metric f^ u = d^^dy^P^AB in spherical co- 
ordinates, this gives a matrix with the only non-zero entries in the upper-left 2x2 matrix, 
and it can be easily reparametrized by using only two Stiickelberg fields. This illustrates 
our point that the reduced massive gravity with the number of scalar fields N less than 
the space-time dimension is equivalent to dRGT theory with a degenerate reference metric 
fftv (or Jab equivalently) . However, the spherically symmetric ansatz given above reduces 
to the degenerate reference metric only in the absence of perturbations. 

3.1 Number of degrees of freedom 

We would now like to estimate the total number of degrees of freedom propagated by the 
full non-linear theory of gravity and two scalar fields. For this we will use the Dirac's 
approach to the Hamiltonian analysis of constrained systems [24, 25]. 

As was already mentioned, due to the fact that the action (3.2) is reparametrization 
invariant and that the scalar fields are coupled to gravity minimally, i.e. only through the 
terms ' d^cj) A d v cj) B , it is legitimate to count the number of degrees of freedom propagated 
by the scalar field action and the Einstein-Hilbert action separately. In such a diffeomor- 
phism invariant theory of gravity and minimally coupled scalar fields, the Hamiltonian 
vanishes on the constraint surface, and both the lapse and the shift enter the Hamiltonian 
linearly. This implies the appearance of in total 2(d + 1) first-class constraints, which can 
be used to reduce the number of gravitational degrees of freedom to d(d+ l)/2 — 2(d + 1). 
The dynamics of the scalar fields then shall be generated by the usual Hamiltonian of the 
scalar field action alone, contained in the Hamiltonian constraint of the full theory. There- 
fore, the scalar field dynamics in such a theory can be considered separately from gravity. 

2 We note that this is not the ansatz usually studied in the context of the spherically symmetric solutions 
of dRGT theory. Instead the common ansatz is <j>° = g(t,r), (p 1 — f(r, t)x l jr. 
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Then naively one would expect that the number of degrees of freedom propagated by any 
dRGT-type massive gravity in (d + l)-dimensional space-time (with d>2) equals to 



where the first term accounts for the degrees of freedom propagated by the massless gravi- 
ton, and the second term is just the number of scalar fields. This naive counting demon- 
strates why, in (3+l)-dimensional space-time, a general non-linear massive gravity theory 
with four Stiickelberg fields propagates six degrees of freedom. It has been demonstrated 
that in the dRGT subclass of massive gravity theories at most five degrees of freedom prop- 
agate due to the special structure of the graviton mass term (2.2) [7, 8]. In Stiickelberg 
language it is clear that, in order for the assertion to be true, the scalar field Lagrangian 
(2.3) has to have a very special structure such that it propagates less degrees of freedom 
than the number of fields. In other words, in the non-linear dRGT massive gravity the four 
Stiickelberg fields do not correspond to four independent degrees of freedom [10] (see also 
[11, 12]). This can be seen from the vanishing of the determinant of the kinetic (Hessian) 
matrix of the scalar field Lagrangian 



Hence the equations of motion of the scalar fields are not independent from each other, and 
there exists (at least) one combination of the equations of motion which gives a constraint 
equation relating the canonical momenta of the scalar fields. As a result, in dRGT massive 
gravity the scalar fields propagate at most iV — 1 = 3 degrees of freedom. 

Our ultimate goal is to find the constraint structure of the scalar field part of the full 
dRGT massive gravity while keeping the space-time metric arbitrary. In this paper we 
start with the case of the reduced massive gravity (3.2) with two scalar fields. For this 
we explicitly calculate the determinant of the kinetic matrix of the theory. Curiously, we 
show that the naive expectation, that also in the case of two scalar fields the determinant 
vanishes and the theory propagates N — 1 = 1 degree of freedom, is not met. Instead we 
find that, in general, the determinant is not equal to zero, and thus there are two dynamical 
degrees of freedom in the scalar field sector. 

3.2 Determinant of the kinetic matrix 

In the case of two scalar fields, the only non-vanishing characteristic polynomials of the 
square root matrix can be explicitly expressed in terms of the TrX and detX as 




(3.6) 



Aab = 




d 2 c, 



(3.7) 





(3.9) 



Then the scalar field action (3.2) reads 




(3.10) 
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Since the 0q term does not affect the dynamics of the scalar fields, in what follows we set 
00 = 0. We also note that in dRGT massive gravity case, where the number of scalar fields 
N coincides with the number of space-time dimensions, the highest order term with (3n is 
usually dropped since it is a total derivative. In the reduced massive gravity, however, the 
term with 0n=2 does contribute to the dynamics of the scalars and, in general, cannot be 
neglected. 

In order to separate the time derivatives of the scalar fields while keeping the space- 
time metric arbitrary, we employ the ADM formalism [26]. In ADM variables for the metric 
components 
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1 

' N 2 



A' ' 

AT2 



N i ij _ N'N' 
N 2 ' N 2 



(3.11) 



the matrix X4 can be expressed as 



X 



B 



D(f> D<f> + S AC ) f BC 



(3.12) 



where D = jj (do — N l di) , and the matrix S AC = 7 JJ di<j) A djcj) C depends only on the spatial 
derivatives of the scalar fields. The canonical momenta conjugated to the scalar fields are 
given by 



1 dC$ 

NdD4> A 



( 



0i 



^(TrZ+2Vdetx) 



1/2 



D<I>a + 



VdetX 



(S Jab - Sab) D0 



b 



+ 



+ 



02 



VdetX 



(Sf A B-S A B)D(j ) l 



(3.13) 



Here S = Tr Sg, and the TrX and detX also depend on the time derivatives Dcf> A as 

TrX = TrX A = S - D<j) A D^> A , 
detX = detX^ = det S - D^ A Dcp B (5 Jab - S A b) , 

where det S = det S A = det / det S AB . The determinant of the kinetic matrix is given by 



det Ab = ~ det g 



det S 
(detX) 2 



0\ 



+ 



TrX+ 2 VdetX 
X A (Tr X + 25) (Sf -S6%)- 25 det 5 

0102 



VdetXX A (3S5% - 25 



+ 



TrX + 2\/detX 



3/2 



VdetX (5 TrX + 4 det 5) 



(25X^ + 5^Tr X) (5f - S5%) - 25 det s] + 0\ det 5 . (3.14) 
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This expression is valid for any choice of the scalar field metric /ab(<^) as long as it does 
not involve the time derivatives of the scalar fields. The determinant depends on temporal 
(contained within the matrix I A ) and spatial derivatives of the scalar fields. In general it 
has a non-zero value which depends on the choice of initial conditions. The only special case 
when the determinant vanishes identically is if we are considering a two-dimensional space- 
time where the matrix S AB = 7 11 9i^3i</> B is a matrix of rank one, and det 5 = 0. This 
case corresponds to the two-dimensional massive gravity and our findings are in agreement 
with the previous work by de Rham et al. [10] . If the det S factor appears also in the 
theory with four scalar fields, then the full dRGT massive gravity in 3 + 1 dimensions also 
has the identically vanishing kinetic matrix, and thus at most five degrees of freedom in 
total. 

We thus conclude that in general the action (3.10) describes two independent dynamic 
fields. However, on the surface det S = in the configuration space, the Lagrangian 
equations of motion are degenerate and determine the second time derivative only for one 
independent combination of fields. In the theory of partial differential equations the solu- 
tions that entirely belong to the detS" = subspace are called singular solutions (cf. [27]). 
In other words, singular solutions of a system of differential equations are the solutions 
which belong to the surface where the number of independent highest time derivatives is 
less than the number of the fields. Such solutions in general are the envelopes of families 
of regular solutions of the system, and at each fixed moment of time coincide with some 
regular solution (or the whole family of regular solutions). It means that the initial condi- 
tions on this surface do not specify a unique solution since there are other solutions of the 
theory which are touching the det S = surface at the initial moment of time. We note 
that this discussion holds only classically. In the full quantum theory the vanishing of the 
determinant of the kinetic matrix signals that the scalar fields are infinitely strongly cou- 
pled, and the quantum effects are crucial for the dynamics of the system near the singular 
surface. 

It is interesting to note that the trivial solution (3.3) with <j> A = x A is on the surface 
det S = 0. However, any perturbations around this solution defined as \ A = 4> x 
will no longer be on the singular surface and will propagate two degrees of freedom. In 
order to understand the dynamics of such field perturbations it is instructive to study the 
behaviour of the system in a close vicinity of the singular surface. Note that for fa ^ the 
condition det S = is also a necessary condition for det ^4 = 0. Therefore in our discussion 
of singular solution we will focus on the singular surface det S = 0. Although in general 
the determinant of the kinetic matrix could vanish also in some other regions of the phase 
space. 

3.3 Hamiltonian analysis away from det S = 

For any initial conditions away from the surface det S = the expression (3.13) for momenta 
is invertible, and the system contains two propagating degrees of freedom. In order to 
qualitatively understand the dynamics of the system in the vicinity of the singular surface 
we fix the scalar fields metric to be flat Jab = Vab and construct the Hamiltonian for the 
limiting cases, when only one of the terms in Lagrangian (3.10) is present. 
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First, we consider the case when /?2 = and /3% = 1. The action with only f5\ term 
present, in the case of four scalar fields, was already studied as a special case of the dRGT 
theory and is named as "minimal non- linear massive gravity". Our Hamiltonian is in 
agreement with the previous results (cf. [12]). The expression (3.13) for the momenta can 
be inverted to give 

1/2 



TrX + 2 v / detX 



/~fu s i dct5\ V VdetX 



^ + —^5^) . (3.15) 



L>0 A 

- ( i -L 

VdetX 

It still does not allow to express the velocities in the terms of momenta completely, but it 
turns out to be enough in order to obtain the Hamiltonian in terms of S AB and ir A . After 
some algebra and with the help of the relation = 1 +7 _1 7r yl 7r^4, the Hamiltonian takes 
the following form: 

1/2 



n = -Ny/j ( S + -f- 1 Tr A 7t B S AB + 2JdetS (1 + -f^ 1 7r A 7r A ) ) +N i d i (i> A K A . (3.16) 



This Hamiltonian has the form % = NT~Lq + N l T~Li, linear in the ADM lapse and shift, 
as it should be in any minimally coupled theory where the scalar fields enter the action 
only through different combinations of gt 11 'd^4> A d u (f> B . When considered together with 
gravity, Ho and Hi simply contribute to the Hamiltonian constraint and to the generators 
of spatial diffeomorphisms respectively. Doing so does not change the dynamics in the 
scalar field sector, and merely reflects the reparametrization invariance of the action. We 
therefore feel free to consider the scalar fields separately from gravity. At last, we note 
also that the scalar fields Hamiltonian T~Lq can be written as a trace of the square root 



matrix H = - Tr J S A {5% + 7- 1 ir c ir B ) , very similar to the structure of the Lagrangian 
In the case when /3i = and /?2 = 1 the velocities can be expressed as 



A Vdetl AB 

D(p = "yodels 5 **■ (3 - 17) 



Using the relation g^-j = 1 H — - d pj 4 g~ B one can obtain the scalar fields Hamiltonian as 



n = - V7 Vdet S + 7- 1 7r A 7T B S AB = -^detJs A (s% + ) ' ( 3 - 18 ) 

The form of this Hamiltonian is also similar to the form of the original Lagrangian, 
£-4> = V~9 det V^j and can be written as a determinant of some square root matrix. 

In order to look at the dynamics we focus on the former case with /3 2 = 0. The 
equations of motion for the scalar fields in Hamiltonian form read 

D<j> A = — ( dGtS = Tr A + S ab tt b ) , (3.19) 

KA-d i {N i ir A )=d i (^- VAB + 7- 1 ir A ir B + ^det 5 (1 + 7- 1 Tr A n A ) S AB d^A, 

(3.20) 
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where the inverse of the spatial derivative matrix is = (Swab — Sab)/ det 5. These 
equations of motion describe the evolution of the scalar fields 4> A and their conjugated 
momenta tta for any initial conditions with det S ^ 0. In general one expects that all the 
regular solutions, i.e. the solutions specified with the initial conditions with det S / 0, are 
tangential to the surface det S = at some point of time. In other words, the singular 
solutions, for which det S = at any time, are the envelopes of the families of regular 
solutions. Choosing the conditions in vicinity of the singular surface and following the 
infinitesimal evolution in time, one can study the phase portrait of the system near the 
singular surface and the way regular solutions are connected to the singular ones. We also 
note that the Hamiltonian (3.16) cannot be used to study the singular solutions themselves 
since it relies on the assumption det S ^ 0. The solutions with det S = shall therefore be 
studied separately. 

3.4 Time evolution in the vicinity of the singular surface 

In order to illustrate the behaviour of regular solutions in the vicinity of the det 5 = 
surface let us choose some initial conditions that are infinitesimally close to the known 
trivial solution cff = t, (f) 1 = x 1 , but have non- vanishing det S and its time derivative. 
Simplest way to write such initial conditions is to provide a small x 2 (or x 3 ) dependence 
to the 0°: 

0°(*o) = t + e x 2 , 4>°(t ) = l + e o x 2 , 4> 1 (t ) = x 1 , <j> 1 (t ) = Q, (3.21) 

where eo, and eo are arbitrary constants, vanishing in the case of the critical solution 
4> A = x A . Since det S = — 6q and ^detS = — 2eo£o, these constants characterize the 
displacement from the singular surface and its time derivative at the initial moment. Using 
the Hamiltonian equation (3.20) for the momenta one can follow the infinitesimal evolution 
of fields (f) A in time. Moreover it happens to be possible to find an exact solution in the 
case of initial conditions (3.21). It can be obtained by promoting the x 2 dependence of the 
initial conditions to be valid at all times. By plugging the ansatz (jP(t) = £(t) +e(t) x 2 into 
the equations of motion (3.20) one obtains two equations for the functions and e(i) 

C(i) = 2^)^||, e(t) = 2^. (3.22) 

The general solution for the fields 4> A (t) is given by 

o (t) = ^-*o) + eO * 2 0i( t )=x 1 . (3.23) 

It happens that this family of solutions never approaches the singular trivial solution (jp = t 
independently of how close are the initial conditions to it, i.e. how small is eo- Instead, 
the solutions (3.23) are asymptotically approaching a different set of singular solutions 
(f)° = £o |^ + £o = const in the limit t — > ±oo. Therefore, for any given constant there is 
a three parameter subfamily of regular solutions that approach it in the t — > ±oo limit. 
Figure 1 illustrates this behaviour for the singular solution <fp = 0. For simplicity we 
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have suppressed the x 2 dependence of the regular solutions, and each line on the figure 1 
corresponds to the one parameter family of solutions, which are different from each other 
by the constant rescaling of e(t). From the solution (3.23) one can see that det S oc h, 
and therefore the det S = surface cannot be reached along the discussed trajectory at 
any finite time. Note that apart from the solutions that start at the finite distance from 
the det S = surface and approach it in the future there exist solutions that start as a 
small perturbations of the singular solution and leave the det 5 = surface. We would also 
like to remark that all the solutions with a given x 2 dependence have a singularity at the 
finite time t = to + Hopefully there are other solutions in this theory that are free of 
singularities. 



4 The singular surface with det S = 

In this section we discuss the most general scalar field configurations which satisfy the 
condition det 5 = and show that the dynamics of the scalar fields in this subspace are 
equivalent to the dynamics of scalar fields in the case of 1 + 1 space-time dimensions. To 
see this we first discuss the 1 + 1 dimensional case separately. 

4.1 1+1 dimensions: massive gravity 

In 1 + 1 dimensions the scalar field Lagrangian (3.10) reduces to 

= PiV^g (Til + 2V det iy /2 , (4.1) 

since the fii term is a total derivative term. The above Lagrangian coincides with the 
dRGT mass term and was previously analyzed in [10, 28]. Here we shall follow a different 
approach of Hamiltonian analysis which enables us to find the gauge symmetry of the scalar 
field action. We show that the scalar field action of non-linear massive gravity propagates 
no degrees of freedom in 1 + 1 dimensions, in agreement with [10, 28]. 
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We first note that in this case the determinant of the matrix is a full square 

detX^ = ^ (^iABE^d^d^A 2 , (4.2) 

where the bared Levi-Civita tensors e^, eab denote the flat space antisymmetric tensors 
defined as e 01 = — e 10 = 1, etc. in every coordinate frame. In this section for simplicity 
we will consider the flat Minkowski scalar field metric Jab = Vab for which the factor 
det/ = — 1. The scalar field action then becomes (up to a constant factor) 



C<h = 2 



N 



1/2 



Drp + + vWA Dij)- - vW_ » ( 4 - 3 ) 



1/2 



where ip± = db <j) , D = (do — N 1 d\)/N, and ijj± = &iip±. In the following we perform 
the full Hamiltonian analysis of this system according to the constraint analysis proposed 
by Dirac and extended by Henneaux et al. [24, 25]. 

4.2 Minkowski background 

We start with the case of a flat Minkowski background metric = rjuv since the gener- 
alization to an arbitrary background is straightforward, as we shall see below. 

4.2.1 Constraint algebra 

In flat space the Lagrangian takes the simple form 



^ = 2^ + +v;v^--^-> ( 4 - 4 ) 

and the conjugated momenta to the fields ip± are 



n + = v , tt_ = V . (4.5) 



It is obvious that the momenta are not independent. Instead, they satisfy the primary 
constraint 

C = 7T+ - — = . (4.6) 

The total Hamiltonian density of the theory is obtained by adding the primary constraint 
to the Hamiltonian as 

= vr-V- - — V>+ + ^oC , (4.7) 

where the Lagrange multiplier uq = Uo(t,x) is an arbitrary function of the space-time 
coordinates. For the analysis of the dynamics of the system we define the equal-time 
Poisson bracket as 



- 13 - 



The time evolution of a functional f(t, x) = f(t, x, ip r (t, x), 7Tj(t, x)) is then given by 



For the consistency of the Hamiltonian equations of motion with the Lagrangian equations 
of motion one has to impose an additional constraint to the system, namely that the 
primary constraint is preserved in time. This in general leads either to secondary (and 
tertiary, . . . ) constraints or determines the arbitrary function uo(t,x) [24]. In our case we 
obtain a secondary constraint 



It is straightforward to check that the time evolution of C\ does not imply any new con- 
straints since 



is a spatial derivative of the secondary constraint itself. Since both constraints mutually 
commute, i.e. {Co(x), Co(x')} = {Co(x), Ci(x')} = {Ci(x), Ci(x')} = 0, and since there are 
no further constraints, we conclude that the constraint algebra is closed and our system 
has two first class constraints. Cq is a primary first-class constraint and C\ is a secondary 
first-class constraint. 

4.2.2 Gauge symmetry 

The existence of first-class constraints indicates that there is a gauge symmetry in our 
theory. The purpose of this section is therefore to identify the gauge symmetries of the 
original Lagrangian (4.4) and find the number of degrees of freedom described by it. 

Since the total Hamiltonian (4.7) contains an arbitrary function of space-time coordi- 
nates uo, a given set of initial conditions for the canonical variables r ip t , 7Tj after some time 
interval will evolve to different values of the canonical variables for different choices of uq. 
Any two such set of values describe the same physical state related by a gauge transforma- 
tion. In order to find the generators of the transformation one considers the evolution of 
a given set of initial data over a finite time interval. This is reached by multiple Poisson 
brackets of the canonical variables and total Hamiltonian, each of them transforming the 
system infinitesimally. Hence after a finite time interval two different sets of canonical vari- 
ables obtained from the same initial data will differ by a gauge transformation generated 
by all first-class constraints. It is therefore why all the first-class constraints should be 
put on the same footing and the Hamiltonian should be extended by adding to it also the 
secondary (and tertiary, . . . ) first-class constraints [24]. This makes the full symmetry of 
the theory manifest. In our case the extended Hamiltonian looks like 




(4.9) 




(4.10) 





1 



■0+ + u C + u\C\ 
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where we have introduced another arbitrary function U\(t, x). Under the transformations 
generated by the two constraints the canonical variables q = {ib l , ^i] transform according 
to the law 

q ^ q + 5q , 5q(x) = lq(x), [ dx' [e (x')Co(x') + ei(a/)d(a/)] ) . (4.13) 



This for the transformations of the canonical fields gives 



5ib + = £ , ^- = ---(eo - 4) , ( 4 - 14 ) 



n 2 _ 

while the conjugated momenta stay unchanged. The corresponding extended first order 
action 

S E = d 2 x n + ip+ + 7T--0- - He (4.15) 

is invariant under the above gauge transformations if also the Lagrange multipliers uq, u\ 
transform. Their transformation laws are not of any need in the present work, therefore we 
shall not give their explicit form and instead refer the reader to [25]. Due to the fact that 
in (4.15) we have introduced an additional arbitrary function u\, the equations of motion 
which can be derived from (4.15) do not coincide with the equations of motion following 
from the action St = / d 2 x (^iip 1 — Ht^J or equivalently from the original action (4.4). 
Moreover, the original Lagrangian is not invariant under the gauge transformations (4.14). 
The reason for this is that the extended Hamiltonian formalism introduces an additional 
redundancy in the description. However, the time evolution of the gauge invariant fields can 
be equally well described by both the total Hamiltonian %t an d the extended Hamiltonian 
H E . 

In order to obtain the symmetry of the original scalar field action, one can rewrite the 
transformations (4.14) by expressing the conjugated momenta according to their definitions 
(4.5) and demand that the action remains unchanged. This leads to the following relation 
between the gauge parameters 

eo = \{e' l -e l ) • (4.16) 

Hence the gauge symmetry of the Lagrangian is 

i J; —ih' 
(s' + s)K ^, (4.17) 

^+^^ + + l -{e' -e) . (4.18) 

Since the above symmetry transformation involves both, the gauge parameter e and its time 
derivative, then the number of degrees of freedom in the theory are reduced by two which 
coincides with the total number of first class constraints [25]. It is so, because the gauge 
parameter and its time derivatives are independent functions in the sense of independent 
initial data which can be chosen arbitrarily at the initial moment of time. 3 Another way 



3 A familiar example where exactly the same approach of counting the degrees of freedom can be applied 
is electrodynamics. There the gauge transformation of the vector field — > A^ + <9 M A also involves both 
the gauge parameter A and its time derivative. The constraint analysis of the theory also shows that there 
is one primary and one secondary first-class constraint removing two out of four degrees of freedom. 
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to see that there are no propagating degrees of freedom is by performing the gauge fixing 
in the extended action (4.15). Since there are two constraints on the momenta and two 
gauge symmetries (4.14) on the canonical fields it is evident that the action is pure gauge 
and propagates no degrees of freedom. The same conclusion could have been drawn also 
from the analysis of the Lagrangian equations of motion. 

4.3 Arbitrary background 

The scalar field action in an arbitrary curved 1 + 1 dimensional space-time can be written 

as 



l/ 2 ■ 

V>+ + a + ip + 4>- ~ a -i>'- 



1/2 



(4.19) 



where we have introduced the notations a± = N\f 7 11 =p N 1 . The conjugated momenta are 
defined as 



7T+ 



'7ii 



ip- — a-tp'_ 
tJj + + a + ip' + 



7T_ 



/7ll 



ip + + a+ip' + 
tp- — a_ip'_ 



(4.20) 



and the Hamiltonian analysis of the system can be carried out in complete analogy to the 
case of Minkowski background. The extended Hamiltonian and the closed set of constraints 
can be found to be 



He = a_7r_7//_ — a + — ip', + uqCq + u\C\ 

IT- 



with Ci 



= k + 



7n 

7T_ 



1 

7T_ 



(4.21) 
(4.22) 



As before the constraints Co and C\ are first class constraints and generate the gauge 
transformations of the canonical variables ifj± \-t ip± + 5tp± with 



5ip+ = eo 



-(■yneo-e'i) ■ 



(4.23) 



By inserting them in the Lagrangian (4.19) one obtains the following condition on the 
gauge variables 



£ 7il(a + + a-) - e (d - a^di - 2a'_) 711 + (d - a_3i - 2a '_) 4 = 0, 



(4.24) 



under which the Lagrangian remains invariant under the transformations (4.23). This 
condition can be rewritten in metric components by using the relations 711 = g 00 det g, 
711(0+ + a-) = 2</ = g, and 



(4.25) 



where the only non-zero components of the Levi-Civita tensor are e 01 = — : 



10 = -iV^g)- 1 . 

Unfortunately, for generic background metric it is impossible to solve (4.24) for eo in local 
form. The gauge transformation is therefore in general non-local. 
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4.4 3+1 dimensions 



For the two scalar fields in 3+1 dimensions the determinant of the matrix of spatial deriva- 
tives S AB = ^di4> A dj4> B reads 



detS AB = 



Jmn 



(4.26) 



Hence the condition det S = translates into requirement that the norm of the cross 
product of the spatial gradients of the scalar fields (jP and (j) 1 vanishes. In other words it 
means that both gradients of the scalar fields have to lie along the same spatial direction 
and thus can be used to parametrize only one spatial direction. Therefore the most general 
scalar field configuration satisfying det S = can be parametrized as 



l> = <f> (t, f(t,x 1 )) 



(4.27) 
(4.28) 



where f(t, x l ) is an arbitrary function of space-time coordinates. 

In order to see that this ansatz for the scalar fields makes the dynamics of the 3 + 1 
dimensional theory equivalent to the dynamics of the 1 + 1 dimensional theory it is useful 
to introduce the short hand notations N = N, N l dif — dof = Nf, difdjf^ = jff . 
In terms of these variables the 3 + 1 dimensional field Z AB takes the form 



rAB 



= g^d^ A d u ^ B = —L(dt- N f d f ) 4> A (dt - N f d f ) ^ + j ff d f cf> A d f (t> 1 



= X AB ^~g^d^ A d^ 



(4.29) 



where the tilded indices take the values fi = 0, /. We recognize the tilded variables 
N, N f , jff as the ADM variables of an effective 1 + 1 dimensional metric g& u . Indeed, for 
the components of the effective metric 



g tt , g of = g u d t f + g u d t f , 



gff = g a d t fd t f + 2g u d t fd t f + g^difdjf 



they satisfy 



9 



oo 



1 



9 



of _ 



(4.30) 



N 2> ~ 9ff=¥f -{w) ■ 

As in the 1 + 1 dimensional case, the determinant detX can be rewritten as a full square 

1 2 



detX = detZ^f 



det/ 



det g [2 



1 



eABe^dfi, 



(4.31) 



with (detg) -1 = —'yff/N 2 , and eab denoting the flat space antisymmetric tensors. 
Hence all the terms in the Lagrangian containing the scalar fields can be rewritten in terms 
of an effective two-dimensional metric g^ u . We would like to emphasize that this rewriting 
is merely cosmetic and has the meaning only as the simplification of notations in the scalar 
field action. 
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In order to simplify the analysis of the equations of motion of the scalars, we rewrite 
the integration measure of the Lagrangian density in another coordinate system {x^}, 
where x° = t, x 1 = f(t,x l ), and x 2 = x 2 (x M ), x 3 = x 3 (x^) some arbitrary non-degenerate 
coordinate transformations. In this case the metric components transform according to the 
usual transformation laws, and the components g 00 , g 01 , g 11 coincide with the components 
of the effective 1 + 1 dimensional metric g^ u , jl = {0, /} given above. Hence the Lagrangian 
of the scalar fields can be rewritten in terms of the metric g^ u as 



dx 2 dx 3 



dtdf \/—q — 



1/2 



1/2 



(4.32) 



where a± = Nyfifff + . The variables N, , 7^, used for notational simplicity only, 
can be expressed in terms of the metric g^ v as in (4.30). By comparing this action with 
(4.19) one sees that the only difference is the volume factor and the prefactor \J—gjN ^ 
■v/7ii, which depends on all four space-time coordinates. Under the assumption that the 
volume spanned by x 2 , x 3 is finite, the Hamiltonian analysis of the scalar field dynamics 
coincides with that in section 4.3. 

We thus conclude that the ansatz for the scalar fields (4.27), (4.28) such that the 
condition det S = is satisfied leads to a theory which is equivalent to the 1 + 1 dimensional 
case and thus propagates no degrees of freedom. In Hamiltonian language, on this subspace 
of the scalar field configurations the theory has two first class constraints. 



5 Conclusions 

Any diffeomorphism invariant formulation of massive gravity inevitably contains a number 
of scalar fields minimally coupled to the dynamical metric field and can be viewed as just 
some particular scalar field theory coupled to general relativity. Therefore we argue that 
the Hamiltonian structure and the counting of degrees of freedom can be done for gravity 
and scalar fields separately. In other words, the absence of the sixth degree of freedom in 
the dRGT non-linear massive gravity [7] can be seen as a feature of the scalar field action, 
and can be studied in the scalar field theory given by the dRGT mass term. 

While the full dRGT scalar action contains the number of fields equal to the space-time 
dimension, in this paper we have focused on the reduced case with two scalar fields, which 
coincides with the full theory only in 1 + 1 dimensions. We have calculated the determinant 
of the kinetic matrix d 2 £ ( f ) /d(f) A d(j) B of the non- linear theory and have found that in d > 1 
dimensions it does not vanish for generic initial conditions. Thus in more than 1 + 1 
dimensions both of the fields are, in general, propagating. However there exists a subspace 
of the configuration space where the Hessian is vanishing. It corresponds to the case where 
the coordinate transformation represented by the scalar fields (f) A (x l ) is singular on any 
two-dimensional space-like surface, or, equivalently, when both of the fields depend only 
on one independent space-like direction. In this case the scalars effectively live on the 1 + 1 
dimensional space-time, and the theory is equivalent to the 1 + 1 dimensional dRGT massive 
gravity, where there is only single spatial direction available. For the latter constrained 
theory the full Hamiltonian analysis reveals two first-class constraints which generate one 
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gauge transformation that leaves the action invariant. Since the transformation involves 
two independent parameters, then after fixing the gauge the theory does not contain any 
degrees of freedom. This is in agreement with the previous findings in the 1 + 1 dimensional 
dRGT massive gravity [10, 28]. For the theory in more than 1 + 1 dimensions the effectively 
1 + 1 dimensional solutions with vanishing Hessian correspond to the so-called singular 
solutions. On such a singular solution at each moment in time there exist infinitely many 
other regular solutions of the theory which are tangential to the singular solution, i.e. 
with coinciding 4> A (x l ) and <j) A (x l ). Therefore, there is no choice of initial conditions that 
uniquely specifies such a solution, and any perturbation in the initial conditions leads to 
the regular solution with two degrees of freedom and non-vanishing Hessian. We note that 
our findings do not allow us to draw conclusions about the behaviour of the dRGT-like 
theories with more than two scalar fields, but the proposed method can be extended to 
include arbitrary number of scalar fields. 
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